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The problem of three interacting Kondo impurities is solved within a numerically exact continuous 
time quantum Monte Carlo scheme. A suppression of the Kondo resonance by interatomic exchange 
interactions for different cluster geometries is investigated. It is shown that a drastic difference 
between the Heisenberg and Ising cases appears for antiferromagnetically coupled adatoms. The 
effects of magnetic frustrations in the adatom trimer are investigated, and possible connections with 
available experimental data are discussed. 
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The electronic structure of adatoms and clusters on 
surfaces constitutes one of the most fascinating subjects 
in condensed matter physics and modern nanotechnology 
The scanning tunneling microscopy (STM) or spec- 
troscopy technique allows the study of atomic structure, 
the electronic energy spectrum, and magnetic properties 
of different surfaces at an atomic scale. In particular, 
STM gives the unique opportunity of directly investi- 
gating an essentially many-body phenomenon, namely 
the Kondo effect |^ Q . Earlier only indirect meth- 
ods such as analysis of temperature and magnetic field 
dependencies of thermodynamic and transport proper- 
ties were available ^sl, ll| . Recently STM studies of small 
transition metal nanoclusters on different surfaces have 
been performed, including Co dimers and Cr trimers 
81 on a Au surface, and Co clusters on carbon nanotubes 
9j. The electron spectrum of these nanosystems, in par- 
ticular the existence of the Kondo resonance, turns out 
to be very sensitive to the geometry of the clusters as 
well as to the type of magnetic adatoms. The later can 
be important for nanotechnological fine tuning of surface 
electronic structure. 

The "quantum-corral" type of STM-experiments pro- 
vides an unique opportunity to investigate in detail 
an interplay between the single-impurity Kondo effect 
and interatomic magnetic interactions in nanoclusters 
which is a key phenomenon in Kondo lattice physics 
m El El El Hi- The interaction between itinerant 
electrons and localized ones leads to the screening of the 
magnetic impurity moment which is the Kondo effect; 
on the other hand, the RKKY exchange interaction be- 
tween localized spins suppresses the Kondo resonance at 
the Fermi level. As a result, a very complicated phase di- 
agram can be obtained with regions described by a strong 
coupling regime, "normal" magnetic behavior with loga- 
rithmic corrections, and non- Fermi-liquid behavior 
Quantum critical points at the bound ary of different 
phases is a subject of special interest 0, There 
exists a general belief that anomalous features of many 
/-electron systems such as heavy-fermion or non-Fermi- 



liquid behavior can be treated in terms of the Kondo 
lattice picture 6, l^L I1(tI|. This is why non-perturbative 
investigations of the basic physical features of few-atom 
magnetic clusters in a metallic medium is of primary in- 
terest. At the same time, due to the extreme complex- 
ity of the problem, theoretical investigations of electronic 
structure for several Kondo centers usually involve some 
uncontrollable approximations, such as a replacement of 
the Heisenberg interatomic exchange interactions by the 
Ising ones jl4j or a variational approach based on a sim- 
ple trial function [l^ . 

In this Letter we present results of a numerically ex- 
act solution of the three Kondo impurity problem within 
the recently developed continuous time quantum Monte 
Carlo (CT-QMC) method [ijI. For the antiferromagnetic 
(AFM) exchange interatomic interaction, in contrast to 
the ferromagnetic (FM) one, the results for the Heisen- 
berg and Ising systems differ essentially. Based on our 
theoretical analysis, the recent paradoxical experimen- 
tal results 8] where the Kondo resonance is observable 
for an isosceles magnetic triangle but not for the perfect 
Cr-trimer or individual Cr adatom will be discussed. 

We start with the system of three impurity correlated 
sites with Hubbard repulsion [/ in a metallic bath and 
with an effective exchange interaction Jij between them, 
a minimal model which however includes all relevant in- 
teractions necessary to describe magnetic nanoclusters on 
a metallic surface. The effective action for such cluster 
in a metallic medium has the following form: 



S = Sq + W, 



(1) 



So 



The last term in the right-hand-side of Eq.Q allows us 
to consider the most important "Kondo lattice" feature, 
that is, the mutual suppression of the Kondo screening 
and intersite exchange interactions . Another fac- 

tor, the coherence of the resonant Kondo scattering, is 



2 



taken into account by the introduction of inter-impurity 
hopping terms Uj to the bath Green function which is 
supposed to be ^^^^ = Q^^^ij ^ Here CJri(icj„) = 
fi + i{u)n + \/^n + 1)72 corresponds to the semicircular 
density of states (DOS) with band-width 2 and are 
inter- impurity hopping integrals. For real adatom clus- 
ters the exchange interactions are mediated by conduc- 
tion electrons (RKKY interactions) which are dependent 
on the specific electronic structure of both adatoms and 
host metal. To simulate this effect we will consider 
as independent parameters which is a common practice 
in the Kondo lattice problem 0, Ull Il4] ; otherwise for 
the half-filled non-degenerate Hubbard model used in our 
calculations the exchange is always antiferromagnetic. In 
the model the geometry of the problem is specified by 
the values of exchange integrals Jij and hopping parame- 
ters Uj . We will concentrate on the case of equilateral 
triangle when Jij — J and = t] to compare with the 
experimental situation in Ref|8| also an isosceles triangle 
will be considered. To check an approximation used in 
Refill we will investigate the case when all spin-fiip ex- 
change terms are ignored and the Heisenberg (SS) form 
of interaction JijSi(r)Sj(T) is transformed into the Ising 
(5,5.) one J,,St{r)SI{r). 

We use the numerically exact CT-QMC method for 
our computer simulations. Unlike the Hirsh-Fye discrete- 
time scheme 20] it docs not involve auxiliary Ising spins, 
but performs a random walk in the space of terms of 
the perturbation expansion for the Green function. One 
of the advantages of this novel approach is the opportu- 
nity it provides to study systems with non-local (in space 
and in time) interactions, which in the usual Hirsh-Fye 
approach would involve a huge increase in the required 
number of auxiliary fields and time slices 0| . 

In brief, numerical simulations make use of the divi- 
sion of the action into to Gaussian and interaction parts. 
The later can be presented in the form J dr'^ Wiji'ji (r), 
where w = Uijiiji (c| ^Cjj — a-^){c[, ^Cj'^i — a^), where the 
a's are c-numbers chosen in a special way to minimize 
the sign problem in the QMC simulation. It is worth 
noting that only the presence of these a' s makes possible 
the simulation of systems with repulsive forces [ig . 

The formal interaction-representation expansion in 
powers of the interaction for the Green function Gij (r, r') 
reads 
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(2) 

Here Z is the partition function and / denotes the set 
Our algorithm performs the random-walk in a 
space of all possible values for k;Ti...Tk;Ii...Ik- The nu- 
merical averaging of Eq. ^ over this random walk gives 
the desired Green function for the interacting system 
with the action (Q. Typical value for k in our calcu- 




FIG. 1: Imaginary part of the Green function at Matsubara 
frequencies for the model ^ with Q~^{iUn) = fi + icJn for SS 
and SzSz interactions. Symbols are CT-QMC data, lines are 
exact diagonalization (ED) results. Parameters of the model: 
(7 = 2, J = 0.2, t = 0, /3 = 8, M = U/2. 



lation is k ^ J dTj2\\wiiT)\\ ^ NPU + N'^PJ, where 
f3 is the inverse temperature and N is number of atoms 
in the cluster. One can see that for the case U ^ NJ 
the exchange interaction indeed does not slow the cal- 
culation down, its complexity is determined by the local 
(Coulomb) interaction. 

In order to check the CT-QMC algorithm for a system 
with complicated Heisenberg interactions we apply this 
method to a simple Hamiltonian analogue of the model 
(Q, i.e. for Q^^{iujn) = M + i^n- We compare our CT- 
QMC approach with the solution obtained using the ex- 
act diagonalization method. Results for the system with 
AFM J (J > 0) are shown in Fig.l for S^S^ and SS in- 
teractions. The estimated errorbar in numerical data is 
10"'^ or less. 

Although the problem of DOS-calculations from G(t) 
data is ill-defined, quite reliable estimations can be made 
at our level of the numerical accuracy. For exam- 
ple, DOS of the above-mentioned Hamiltonian models 
is at most contributed by the four (5-peaks located at 
±U/2, ±{U/2 + J) for Ising and ±C//2, ±(U/2 + 2 J) for 
Heisenberg interaction. Fit of the numerical data with 
several (5-peaks indeed resolves their positions (with a 5% 
errorbar) and relative heights (with a 20% errorbar) . Fur- 
ther, we study models with a continuous DOS, and stan- 
dard maximum-entropy analytical continuation method 
(21j is used to recover DOS. Normally it resolves the DOS 
features from imaginary time QMC data with a similar 
accuracy. 

Let us discuss correlated adatom trimer in the metallic 
bath depending on type of the effective exchange inter- 
action {S^Sz or SS) for AFM and FM cases. First we 
show that SS type of interaction suppresses the reso- 
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FIG. 2: Imaginary part of the Green functions at Matsub- 
ara frequencies for the correlated adatom equilateral triangle 
in the metallic bath for AFM (upper figure) and FM (lower 
figure) types of effective exchange interaction. Parameters: 
U = 2,J = ±0.2, t = 0,/3 = 16,Ai = U/2. There are three 
dependencies on each picture for SS, SzS^ and J = (which 
corresponds to single atom in the metallic bath) types of in- 
teraction. The insets show DOS. 



nance in AFM case. We study the equilateral triangle 
at half-filling. The Green functions at the Matsubara 
frequencies obtained by the CT-QMC technique and cor- 
responding DOS are presented in Fig. 2. The case of 
J = corresponds to the single Kondo impurity in the 
bath. One can see that there is no essential difference 
between S^Sz and SS types of interaction in FM but 
for the AFM one this difference is very important. The 
AFM SS interaction leads to pronounced suppression of 
the Kondo resonance at Fermi level for physically rele- 
vant values of J (Fig. 2). On the other hand, in the FM 
case we do not observe any essential difference between 
the SzSz and SS types of the interaction for a wide range 
of the model parameters. 

To explain these results we need to calculate the spec- 
tral density D{uj) of the on-site spin- flip operators Sf^ for 



an AFM equilateral triangle I?ss(w) = \5{w) + \5{ijJ — 
2>J), Ds^sA'^) = i'^('^) + |'^('^~2J). In both cases there 
is a part of the spectral density with zero frequency due 
to degeneracy of the ground state and spin-flip transi- 
tions between its components, but for the SS case this 
part is twice as large. For the case J ^ Tk [Tk is one- 
site Kondo temperature) only this "soft" component of 
the spectral density will lead to Kondo screening which 
means that the suppression of the Kondo effect is twice 
more efficient for the SS case than for the SzSz one. Since 
the trimer as a whole has a degenerate ground state there 
is still a strong-coupling regime and an effective Kondo 
temperature Tj*. which however is much smaller than Tk. 
Scaling considerations similar to one proposed in Refllll 
gives an estimation of ~ T^l"^ jj^l"^ and ~ T\l P 
for the SS and SzSz case respectively. We assume that 
this quantity is too small to be visible in our simulations 
(as well as in the experimental data Q) so what is ob- 
served corresponds to the one-site Kondo resonance at 
the condition Tk > J ■ Similar estimation for the case 
of FM interactions shows that there is no difference be- 
tween the SzSz and SS model there and ~ T^kI ^'^ 
in the AFM SzSz case. 

In order to describe the experimental situation we 
changed the geometry of the adatom trimer. An observa- 
tion of the Kondo resonance reconstruction was reported 
for one isosceles geometry of three Cr atoms on a gold 
surface j3] ■ Thus we study the isosceles triangles for AFM 
and FM types of effective exchange interaction. We have 
chosen the following parameters of Jij to imitate the ex- 
perimental system: J23 = J, J12 = Jji^Jxz = J/3. The 
computational results are presented in Fig. 3, where one 
can see the reconstruction of resonance in AFM and FM 
cases in accordance with experimental data. Note that 
the Kondo resonance appears only for the more weakly 
bonded adatom in AFM case. 

The observed picture can be drastically changed by in- 
troducing a non-zero value of the hopping parameter tij . 
As noted above the trimer ground state is degenerate at 
tij = 0, however this degeneracy is lifted for tij ^ 0. If 
one of the obtained states lies below the Fermi level, DOS 
can be changed drastically and at certain parameters it 
leads to the appearance of a resonance on the Fermi level 
(see Fig. 4). In Fig. 4, DOS is shown for various values of 
filling in the system and at nonzero tij. It is necessary 
to point out that the introduction of the parameter tij, 
the variation of the filling in the system (/x), and different 
geometries ( Jy ) can lead to various results in DOS. This 
can be used as one of the possible explanations of the ex- 
perimental data regarding the Cr trimer on a Au surface 
which look initially appear counter-intuitive: the Kondo 
temperature for the trimer is much larger than for the 
single site, despite the suppression of the Kondo effect 
by Jij . One can assume that this is a consequence of the 
change of the number of d-electrons in ground-state con- 
figuration for the Cr atom in the trimer in comparison 
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FIG. 3: DOS for equilateral triangle (ET) and isosceles trian- 
gle (IT) geometries with AFM (upper figure) and FM (lower 
figure) types of effective exchange interaction. Parameters are 
the same as in Fig. 2. Values of the effective exchange inte- 
grals for IT are as follows: J23 = J, J12 = J/3, J13 = J/3. 
There are two dependencies in case of IT: one for adatom 1 
and another for equivalent adatoms 2 and 3. All adatoms are 
equivalent in the case of ET (one dependence). 



with the isolated one. In order to describe specific experi- 
mental results quantitatively the method developed here 
should be combined with first-principle calculations in 
the spirit of recently developed "local density approxima- 
tion plus dynamical mean-field theory" (LDA + DMFT) 
scheme |22l| . 

In conclusion, we have shown that the electronic struc- 
ture of a correlated adatom trimer on a metallic surface 
drastically depends on the symmetry of magnetic inter- 
actions. The effective exchange interaction of SS type 
leads to more efficient suppression of the Kondo reso- 
nance in the AFM case than in the case of SzSz inter- 
actions. The experimental STM data can be repro- 
duced qualitatively well by variation of the geometry of 
the problem, hopping integrals and electronic filling for 
magnetic nanosystems. 

The work was supported by POM project N0703M, 
NWO project 047.016.005 and "Dynasty" foundation. 
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FIG. 4: DOS for equilateral triangle with an AFM type of 
effective exchange interaction at various values of filling in 
the system. Parameters: U = 2,J = 0.2, t — 0.3, /3 — 16. 
Corresponding numbers of particles are n — 0.8, n — 1 (half- 
filled case) and n — 1.3. 
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